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We argue that at moderately large momentum transfer —t < lOGeV 2 , hadronic form factors 
and wide-angle Compton scattering amplitudes are dominated by mechanism corresponding to 
overlap of soft wave functions. We show that the soft contribution in both cases can be de- 
scribed in terms of the same universal nonforward parton densities (ND's) J-{x;t), which are the 
simplest hybrids of the usual parton densities and hadronic form factors. We propose a sim- 
ple model for ND's possessing required reduction properties. Our model easily reproduces the 
observed magnitude and the dipole t-dependence of the proton form factor Ff(t) in the region 
1 GeV 2 < — t < lOGeV 2 . Our results for the wide-angle Compton scattering cross section follow 
the angular dependence of existing data and are rather close to the data in magnitude. 
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I. INTRODUCTION 



The Compton scattering in its various versions provides a unique tool for studying hadronic structure. The 
Compton amplitude probes the hadrons through a coupling of two electromagnetic currents and in this aspect it 
can be considered as a generalization of hadronic form factors. In QCD, the photons interact with the quarks of 
a hadron through a vertex which, in the lowest approximation, has a pointlike structure. However, in the soft 
regime, strong interactions produce large corrections uncalculable within the perturbative QCD framework. To 
take advantage of the basic pointlike structure of the photon-quark coupling and the asymptotic freedom feature of 
QCD, one should choose a specific kinematics in which the behavior of the relevant amplitude is dominated by short 
(or, being more precise, lightlike) distances. The general feature of all such types of kinematics is the presence of 
a large momentum transfer. For Compton amplitudes, there are several situations when large momentum transfer 
induces dominance of configurations involving lightlike distances: 

i) both photons are far off-shell and have equal spacelike virtuality: virtual forward Compton amplitude, its 
imaginary part determines structure functions of deep inelastic scattering (DIS); 

ii) initial photon is highly virtual, the final one is real and the momentum transfer to the hadron is small: deeply 
virtual Compton scattering (DVCS) amplitude; 

Hi) both photons are real but the momentum transfer is large: wide-angle Compton scattering (WACS) amplitude, 
the study of which is the ultimate goal of the present paper. 

Our main statement is that, at accessible momentum transfers \t\ < 10 GeV 2 , the WACS amplitude is dominated 
by handbag diagrams, just like in DIS and DVCS. In the most general case, the nonperturbative part of the handbag 
contribution is described by nonforward double distributions (DD's) F(x, y; t), G(x, y; t), etc., which can be related 
to the usual parton densities f(x), A/(x) and nucleon form factors like F±(t), Ga(£)- Among the arguments of 
DD's, x is the fraction of the initial hadron momentum carried by the active parton and y is the fraction of the 
momentum transfer r. The description of WACS amplitude simplifies when one can neglect the y-dependence 
of the hard part and integrate out the y-dependence of the double distributions. In that case, the long-distance 
dynamics is described by nonforward parton densities (ND's) J-(x; t), G(x; t), etc. The latter can be interpreted as 
the usual parton densities f(x) supplemented by a form factor type i-dependence. We propose a simple model for 
the relevant ND's which both satisfies the relation between T{x; t) and usual parton densities f(x) and produces a 
good description of the Fi(t) form factor up to t ~ —10 GeV 2 . We use this model to calculate the WACS amplitude 
and obtain results which are rather close to existing data. 

II. VIRTUAL COMPTON AMPLITUDES 

The forward virtual Compton amplitude whose imaginary part gives structure functions of deep inelastic scat- 
tering (see, e.g., 0) is the classic example of a light cone dominated Compton amplitude. In this case, the "final" 
photon has momentum q' — q coinciding with that of the initial one. The momenta p,p' of the initial and final 
hadrons also coincide. The total cm energy of the photon-hadron system s — {p + q) 2 should be above resonance 
region, and the Bjorken ratio XBj — Q 2 /2{pq) is finite. The light cone dominance is secured by high virtuality of 
the photons: — q 2 = Q 2 > 1 GeV 2 . In the large-Q 2 limit, the leading contribution in the lowest a s order is given 
by handbag diagrams in which the perturbatively calculable hard quark propagator is convoluted with parton dis- 
tribution functions f a (%) (a = u,d, s, . . .) which describe/parametrize nonperturbative information about hadronic 
structure. 

The condition that both photons are highly virtual may be relaxed by taking a real photon in the final state. 
Keeping the momentum transfer t = (p — p') 2 to the hadron as small as possible, one arrives at kinematics of the 
deeply virtual Compton scattering (DVCS) the importance of which was recently emphasized by X. Ji Q (see also 
||). Having large virtuality Q 2 of the initial photon is sufficient to guarantee that in the Bjorken limit the leading 
power contributions in 1/Q 2 are generated by the strongest light cone singularities 0-|?J , with the handbag diagrams 
being the starting point of the a s expansion. The most important contribution to the DVCS amplitude is given 
by a convolution of a hard quark propagator and a nonperturbative function describing long-distance dynamics, 
which in the most general case is given by nonforward double distributions (DD's) F(x, y; t), G(x, y; i), . . . 
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FIG. 1. a) General Compton amplitude; b) s-channel handbag diagram; c) it-channel handbag diagram. 

The DD's arc rather complicated functions. They specify the fractions xp and yr of the initial hadron momentum 
p and the momentum transfer r = p — p' carried by the active parton: k ~ xp + yr. The DD's vanish outside the 
triangle region < x + y < 1 j3fi. In addition to x and y, they also depend on the invariant momentum transfer 
t = (p' — p) 2 . In some limiting cases, the double distributions reduce to simpler and already known functions. For 
r = 0, the matrix elements coincide with the forward ones defining the usual parton densities. This results in the 
following "reduction relations" 



F a (x,y;t = 0)dy = f a (x) 



(1) 



Integrating properly weighted sums of quark and antiquark DD's over x one obtains the Dirac form factor: 

Ve Q /dx/ [F a (x,y;t)-F a (x,y;t)] dy = F 1 (t), (2) 
Jo Jo 

where e a is the electric charge of the "a" quark. Just like for form factors, one should take into account extra 
double distributions K a (x, y; t) corresponding to a hadron helicity flip in the nonforward matrix element These 
distributions are related to the Pauli form factor F 2 (t): one should just substitute F a a by K a a and Fi by F 2 in 
Eq.(||). A common element of these reduction formulas is an integration over y. Hence, it is convenient to introduce 
intermediate functions 



F a (x;t)= F a (x,y;t)dy ; IC a (x;t) = / K a (x,y;t) dy . 

Jo Jo 



They satisfy the reduction formulas 



T a (x- 1 = 0) = f a (x) ; V e a [ [T a {x; t) - T a (x- 1)] dx = F^t) 

Jo 



(3) 



(4) 



V e Q f [K a (x-t)~ /C s (x; tj] dx = F 2 (t) , (5) 

which show that these functions are the simplest hybrids of the usual parton densities and form factors. For this 
reason, we call them nonforward parton densities (ND's). Note that the t = limit of the "magnetic" ND's exists: 
IC a (x;t = 0) = k a (x). In particular, the integral 

y2 e a [k a (x) ~ k a {x)} dx = K p (6) 

gives the anomalous magnetic moment of the proton. The knowledge of the x- moment of fc a (x)'s is needed to 
determine the contribution of the quark orbital angular momentum to the proton spin 0. Since the if- type DD's 
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are always accompanied by the = p^—p'^ factor, they are invisible in deep inelastic scattering and other inclusive 
processes related to strictly forward r — matrix elements. 

There are also parton-helicity sensitive double distributions G a (x,y;t) and P a (x,y;t). The first one reduces 
to the usual spin-dependent densities Af a (x) in the r — limit and gives the axial form factor Fa(£) after the 
x, ^-integration. The second one involves hadron helicity flip and is related to the pseudoscalar form factor Fp(t). 

In the DVCS kinematics, \t\ is assumed to be small compared to Q 2 , and for this reason the t- and m^-dependence 
of the short-distance amplitude in refs. was neglected^ This is equivalent to approximating the active parton 
momentum k by its plus component alone: k — > xp + + yr + . Treating £ = r + /p + as an external parameter and 
using the total fraction X = x + £y as an independent variable, one arrives at an alternative description of the 
DVCS scaling limit in terms of the nonforward parton distributions 0,1) (NFPD's)fl J 7 ^ (X;t). They are related 
to double distributions by 

r mm{x/t,x/0 

T a c a {X ] t)= F a ' a (X-( yi y;t)dy. (7) 

Jo 

In a similar way, one can incorporate the relevant double distributions to define also "magnetic" /C"' Q '(X;t) and 

parton-helicity sensitive nonforward distributions Q^' a (X;t) and V"' a (X;t) p, p^flq| . In addition to the usual 
parton momentum fraction variable X and the invariant momentum transfer t, the NFPD's also depend on the 
skewedness parameter £ = r + /p + specifying the longitudinal momentum asymmetry of the nonforward matrix 
element. This asymmetry appears because it is impossible to convert a highly virtual initial photon into a real final 
photon without a longitudinal momentum transfer. In general, one can use different pairs of vectors to specifiy the 
longitudinal direction: (p, q), (p, q') or (P, q) with P = (j>+p')/2, etc., resulting in different ^-dependent expressions 
for £. However, in the (formal) scaling limit t — > 0, m 2 — » all these expressions for the skewedness parameter £ 
coincide with the Bjorken ratio xsj — Q 2 /2(pq) [||||]- 



III. MODELING ND'S 



Our final goal in the present paper is to get an estimate of the handbag contributions for the large-i real Compton 
scattering. Since the initial photon in that case is also real: Q 2 = (and hence xsj = 0), it is natural to expect 
that the nonpcrturbative functions which appear in WACS correspond to the £ = limit of the nonforward parton 
distributions^] !F^{x; t). It is easy to see from Eqs.(||),(0) that in this limit the NFPD's reduce to the nonforward 
parton densities T a (x;t) introduced above: 

F a Q=0 {x-t)=F a {x-t). (8) 

Note that ND's depend on "only two" variables x and t, with this dependence constrained by reduction formulas 
(H) )(H)- Furthermore, it is possible to give an interpretation of nonforward densities in terms of the light-cone wave 
functions. 

Consider for simplicity a two-body bound state whose lowest Fock component is described by a light cone wave 
function ^(x, k±_). Choosing a frame where the momentum transfer r is purely transverse r — r±, we can write 
the two-body contribution into the form factor as (T(| 

FW(t)=J dx J 9*{x,k ± +xr±)^{.x,k±)^, (9) 



t One should not think that such a dependence is necessarily a higher twist effect: the lowest twist contribution has a 
calculable dependence on t and rrip analogous to the Nachtmann-Georgi-Politzer 0(rrip/Q 2 ) target mass corrections in DIS 

SB- _ m 

'"Other terminology: "off-forward" H, "non-diagonal" |ll| and "off-diagonal" (1211131 is also used in the literature. Off- 
forward parton distributions introduced by X. Ji in his pioneering papers on DVCSHpJ are equivalent though not identical 
to the nonforward ones, while "non-diagonal" and "off-diagonal" distributions essentially coincide with NFPD's, see |Q| for 
details. 

^Provided that one can neglect the t-dependence of the hard part, see a footnote above and discussion in Sec. VI. 
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a) b) 

FIG. 2. a) Structure of the effective two-body contribution to form factor in the light cone formalism, b) Form factor as 
an x-integral of nonforward parton densities. 

where x = 1 — x. Comparing this expression with the reduction formula (^), we conclude that 

^ tb \x,t) = J y*( x ,k x + xrx_)*(x,k ± )^ (10) 

is the two-body contribution into the nonforward parton density !F(x,t). Assuming a Gaussian dependence on the 
transverse momentum k± (cf . JToj ] ) 

V(x,k±.) = $(a;)e- fc i/ 2a;2A2 , (11) 

we get 

jr(t&) t ) = j(tb) ( x -) e st/Ax\* : ^ 

where 

f&\ x ) = ^&{x)=J**>(x t t = ti) (13) 

l07T^ 

is the two-body part of the relevant parton density. Within the light-cone approach, to get the total result for 
either usual f(x) or nonforward parton densities !F(x,£), one should add the contributions due to higher Fock 
components. By no means these contributions are small, e.g., the valence du contribution into the normalization 
of the 7r + form factor at t — is less than 25% |ll| . In the absence of a formalism providing explicit expressions 
for an infinite tower of light-cone wave functions we choose to treat Eq. ([l2|) as a guide for fixing interplay between 
the t and x dependences of ND's and propose to model them by 

TTXXA J 

The functions f a {x) here are the usual parton densities assumed to be taken from existing parametrizations like 
GRV, MRS, CTEQ, etc. In the t = limit (recall that t is negative) this model, by construction, satisfies the 
first of reduction formulas (||). Within the Gaussian ansatz ([l4]), the basic scale A specifies the average transverse 
momentum carried by the quarks. In particular, for valence quarks 

(kl) a = ^ ^ xxr a al (x)dx, (15) 

ly a Jo 

where N u = 2, Nd = 1 are the numbers of the valence et-quarks in the proton. 

To fix the magnitude of A, we use the second reduction formula in (^) relating J- a (x, t)'s to the F i(t) form factor. 
To this end, we take the following simple expressions for the valence distributions 

fl a \x) = 1.89aT°- 4 (l - a;) 3 - 5 (l + 6x) , (16) 
f%*{x) = 0.54x-°- 6 (l - x) 4 - 2 {l + 8x) . (17) 
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FIG. 3. Ratio Ff(t)/D(t) of the Ff{t) form factor to the dipole fit D(t) = 1/(1 - t/0.71 GeV 2 ) 2 . Curve is based on 
Eqs.(16)-(18) with A 2 = 0.7 GeV 2 . Experimental data are taken from ref.[18]. 



They closely reproduce the relevant curves given by the GRV parametrization J17| at a low normalization point 
Q 2 ~ 1 GeV 2 . The best agreement between our model 

F? ft (t) = [ [e u f: al (x) + e d r d al {x)] e st ' ix>? dx (18) 
Jo 

and experimental data |Q in the moderately large t region 1 GeV 2 < \t\ < 10 GeV 2 is reached for A 2 = 0.7 GeV 2 
(see Fig.||). This value gives a reasonable magnitude 

(k 2 ± ) u = (290MeV) 2 , (k 2 ± ) d = (250MeV) 2 (19) 

for the average transverse momentum of the valence u and d quarks in the proton. 

Similarly, building a model for the parton helicity sensitive ND's Q a (x,t) one can take their t = shape from 
existing parametrizations for spin-dependent parton distributions Af a (x) and then fix the relevant A parameter by 
fitting the Ga{£) form factor. The case of hadron spin-flip distributions K. a (x,t) and V a (x;t) is more complicated 
since the distributions k a (x), p a {x) are unknown. 

At t = 0, our model by construction gives a correct normalization Ff(t = 0) = 1 for the form factor. However, 
if one would try to find the derivative [d/ 'dt)Ff(t) at t = by expanding the exponential exp[xt/x\ 2 ] into the 
Taylor series under the integral (|l8|), one would get a divergent expression. An analogous problem is well known 
in applications of QCD sum rules to form factors at small t ]l9|-|22]] . The divergence is related to the long-distance 
propagation of massless quarks in the i-channel. Formally, this is revealed by singularities starting at t = 0. 
However, Ff (t) should not have singularities for timelike t up to 4m 2 , with the p-meson peak at t = m 2 ~ 0.6 GeV 2 
being the most prominent feature of the t-channel spectrum. Technically, the singularities of the original expression 
are singled out into the bifocal correlators |23j which are substituted by their realistic version with correct spectral 
properties (usually the simplest model with p and p 1 terms is used) . An important point is that such a modification 
is needed only when one calculates form factors in the small-t region: for — t > 1 GeV 2 , the correction terms should 
vanish faster than any power of 1/t pjj . In our case, the maximum deviation of the curve for Ff(t) given by 
Eq.(|l8|) from the experimental data in the small-t region -t < 1 GeV 2 is 15%. Hence, if one is willing to tolerate 
such an inaccuracy, one can use our model starting with t — 0. 

IV. SOFT VS HARD CONTRIBUTIONS TO FORM FACTORS 

Our curve is within 5% from the data points Jlq] for 1 GeV 2 < — t < 6 GeV 2 and does not deviate from them 
by more than 10% up to 9 GeV 2 . Modeling the t-dependence by a more complicated formula (e.g., assuming a 



slower decrease at large t, and/or choosing different A's for u and d quarks and/or splitting ND's into several 
components with different A's, etc.) or changing the shape of parton densities f a (x) one can improve the quality 
of the fit and extend agreement with the data to higher t. Such a fine-tuning is not our goal here. We just want to 
emphasize that a reasonable description of the F\{t) data in a wide region 1 GeV 2 < |i| < 10 GeV 2 was obtained 
by fixing just a single parameter A reflecting the proton size. Moreover, we could fix A from the requirement that 
(k\) ~ (300 MeV) 2 and present our curve for Fi(t) as a successful prediction of the model. We interpret this 
success as an evidence that the model correctly catches the gross features of the underlying physics. 

Since our model implies a Gaussian dependence on the transverse momentum, it includes only what is usually 
referred to as an overlap of soft wave functions. It completely neglects effects due to hard pQCD gluon exchanges 
generating the power-law 0((a s /ir) 2 /t 2 ) tail of the nonforward densities at large t. It is worth pointing out here 
that though we take nonforward densities J- a (x,t) with an exponential dependence on i, the F±(t) form factor in 
our model has a power-law asymptotics F-f oft (t) ~ (— 4A 2 /t)" +1 dictated by the (1 — x) n behavior of the parton 
densities for x close to 1. This connection arises because the integral ( |l8| ) over x is dominated at large t by the region 
x r~j AX 2 /\t\. In other words, the large-i behavior of -Fi(t) in our model is governed by the Feynman mechanism 
One should realize, however, that the relevant scale 4A 2 = 2.8 GeV 2 is rather large. For this reason, when 
\t\ < 10 GeV 2 , it is premature to rely on asymptotic estimates for the soft contribution. Indeed, with n = 3.5, 
the asymptotic estimate is F-f oft (i) ~ t~ 4 ' 5 , in an apparent contradiction with the ability of our curve to follow the 
dipole behavior. The resolution of this paradox is very simple: the maxima of nonforward distributions J- a (x,t) 
for \t\ < 10 GeV 2 are at rather low x- values x < 0.5. Hence, the x-integrals producing Ff oft (t) are not dominated 
by the x ~ 1 region yet and the asymptotic estimates are not applicable: the functional dependence of Fj 3oft (i) in 
our model is much more complicated than a simple power of 1/t. 

The fact that our model closely reproduces the experimentally observed dipole-like behavior of the proton form 
factor is a clear demonstration that such a behavior may have nothing to do with the quark counting rules 
Ff(t) ~ 1/t 2 |24|,|25| valid for the asymptotic behavior of the hard gluon exchange contributions. Our explanation 
of the observed magnitude and the ^-dependence of F%(t) by a purely soft contribution is in strong contrast with 
that of the hard pQCD approach to this problem. Of course, there is no doubt that in the formal asymptotic 
limit \t\ — > oo, the dominant contribution to the Fi(t) form factor in QCD is given by diagrams involving two 
hard gluon exchanges, with nonperturbative dynamics described by distribution amplitudes (DA's) tp p (xi, X2, £3), 
(p p (yi, 2/2; 2/3) of the initial and final protons p^j2^ |. However, attempting to describe the data at accessible t by 
hard contributions only, one is forced to make several unrealistic assumptions. 

The crucial element is the use of humpy DA's similar to those proposed by Chernyak and I. Zhitnitsky |28 27 
(CZ). The usual claim is that these DA's are backed by QCD sum rule calculations of their lowest moments. 



However, as we argued in ref. [|29 
because of poor convergence of t 



a straightforward version of the QCD sum rule approach in this case is unreliable 
;he underlying operator product expansion (OPE). In the analysis of the QCD sum 
rules for the moments of the pion distribution amplitude performed in refs. p9|j30[ |, the contribution of exploding 
higher terms of the OPE (neglected in the CZ approach) was modeled by nonlocal condensates. The resulting 
QCD sum rule produces the pion DA close to the asymptotic one. The statement that the pion DA is close to its 
asymptotic form even at a low normalization point is also supported by a lattice calculation of the second moment 
of the pion DA f3l| , by QCD sum rule estimate of the magnitude of ^(x) at the middlepoint x = 1/2 fl32| , by 
the analysis of QCD sum rules for the nondiagonal correlator |33|,|34| , by calculation of the pion DA in the chiral 
soliton model |55| and by a direct QCD sum rule calculation of the large-Q 2 behavior of the 7*777° form factor 
fjfjj. Furthermore, within the light-cone QCD sum rule approach one can relate the pion DA to the pion parton 
densities \p7\ known experimentally. According to the analysis performed in |38| , existing data favor the asymptotic 



shape. Finally, the humpy pion DA advocated in |39|27j] is now ruled out by recent experimental data |4Cj on the 
7*77T° form factor. The data are fully consistent with the next-to-leading pQCD prediction calculated using the 
asymptotic DA 

Since the structure of OPE in the pion and nucleon cases is very similar, we see no reason to expect a significant 
deviation of the nucleon DA from its asymptotic form. In particular, an evidence against humpy nucleon DA's is 
provided by a lattice calculation Q which does not indicate any significant asymmetry. One may argue that the 
proton DA must be asymmetric to reflect the fact that the it-quarks carry on average a larger fraction of the proton 
momentum than the (i-quarks. As shown in ref. to accomodate this observation one needs only a moderate 
shift of the DA maximum from the center point x\ = X2 = £3 = 1/3. Such a shift does not produce a drastic 
enhancement of the hard contribution provided by the humpy DA's. However, with the asymptotic DA, the leading 
twist hard contribution completely fails to describe the data: it gives zero for the proton magnetic form factor 
and a wrong-sign (positive) contribution for the neutron magnetic form factor, with the absolute magnitude of the 
latter being two orders of magnitude below the data |Q. 

Furthermore, as emphasized in refs. the whole strategy of getting enhancements from the humpy DA's 
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is based on an implicit assumption that one may use the perturbative expressions S(k) ~ k/k 2 , D(k) ~ 1/k 2 for 
quark and gluon propagators up to very small virtualities k 2 < (300MeV) 2 . It is worth recalling now why CZ-type 
DA's give an enhanced contribution. Since quarks in the proton carry only a fraction of the proton momentum, the 
characteristic virtualities ~ XiVjt of "hard" quarks and gluons inside the short-distance subprocess are smaller than 
the total momentum transfer t. For a symmetric distribution, one would expect that (xi) ~ 1/3. With the humpy 
DA's, the average Xi for one of the u-quarks is close to 1, and the dominant contribution comes from configurations 
in which this quark is active. Then fractions xi related to passive quarks are rather small. It is precisely the small 
magnitude of the ~ XiXjjt denominators of quark and gluon propagators which produces the enhancement in the 
case of the CZ-type DA's. Hence, to get large hard contributions, it is absolutely necessary to assume that the 
perturbative expressions S{k) ~ fc/fc 2 , D(k) ~ 1/k 2 for quark and gluon propagators may be trusted up to very 
small virtualities. 

An instructive illustration of possible modifications due to finite size or transverse momentum effects is given by 
the light-cone calculation of the 7*777° amplitude |]l6| , |42f in which hard propagator of a massless quark is convoluted 
with the two-body wave function ^(x,k±). Assuming the Gaussian dependence 1 $>(x,k±) ~ exp[— k 2 L /2xxa] on 
transverse momentum, one can easily calculate the k± integral to see that the pQCD propagator factor 1/xQ 2 
is substituted by the combination (1 — exp[— xQ 2 / 2xa\) / xQ 2 which monotonically tends to a finite limit l/2er as 
x — > 0. Hence, the effective virtuality is always larger than 2a. The suppression of propagators at low virtualities 
has a simple explanation: propagation of quarks and gluons in the transverse direction is restricted by the finite 
size of the hadron. Numerically, 2a « 1.35 GeV 2 in that case. However, even a milder modification of the "hard" 
propagators by effective quark and gluon masses 1/fc 2 -> l/(k 2 - M 2 ) with M 2 ~ 0.1 GeV 2 or model inclusion 
of transverse momentum effects strongly reduces the magnitude of hard contributions [ ft9| , especially when the 
CZ type DA's are used. For these reasons, a scenario with humpy DA's and bare ~ 1/xiyjt propagators (which 
amounts to ignoring finite-size effects) considerably overestimates the size of hard contributions. 

The relative smallness of hard contributions can be easily understood within the QCD sum rule context. The 
soft contribution is dual to the lowest-order diagram while the gluon exchange terms appear in diagrams having a 
higher order in a s which results in the usual a s /Tr ~ 1/10 suppression factor per each extra loop. In particular, 
the a s /-K suppression factor is clearly visible in the expression for the hard contribution to the pion form factor 

^r d (Q 2 )w = ^=2(v)$- (20) 

Here, the combination so = 47r 2 / 2 w 0.67 GeV 2 ~ m 2 is what is usually called the "typical hadronic scale" in the 
case of the pion. At asymptotically high Q 2 , the 0(a s /ir) suppression of the hard terms is more than compensated 
by their slower decrease with Q 2 . However, such a compensation does not occur in the subasymptotic region 
where the soft contributions, as we have seen, may have the same effective power behavior as that predicted by the 
asymptotic quark counting rules for the hard contributions. In ref. |54| , both the soft contribution and the 0(a s ) 
corrections for the pion form factor were calculated together within a QCD sum rule inspired approach. The ratio 
of the 0(a s ) terms to the soft contribution was shown to be in full agreement with the expectation based on the 
a s /ir per loop suppression. 



V. COMPTON SCATTERING AMPLITUDE AT LARGE MOMENTUM TRANSFER 

With both photons real, it is not sufficient to have large photon energy to ensure short-distance dominance: 
large-s, small-i region is strongly affected by Regge contributions. Hence, having large \t\ > 1 GeV 2 is a necessary 
condition for revealing short-distance dynamics. 

Consider the Compton scattering amplitude for large values of s, u and t. According to a general rule (see, e.g., 
Q and references therein), to find possible mechanisms generating power-law contributions in the asymptotic limit 
s ~ —u ~ —t ~ Q 2 (Q here is just a characteristic scale), we should look for subgraphs whose contraction into 
point or removal from the diagram kills its dependence on large variables. Contracted subgraphs correspond to 
short-distance (SD) or hard regime while the removed ones to the infrared (IR) or soft regime. Some possibilities 
are shown in FigjiJ 

The power counting estimate for each SD-subgraph is given by 

Ah(Q) ~ Q 4 - N -^ , (21) 

where N is the number of the external photon lines of the hard subgraph H and t$ is the twist of its zth external 
parton line (t = 1 for quarks and physical gluons and t — for longitudinal gluons). 
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a) b) c) d) 

FIG. 4. Some configurations responsible for power-law asymptotic contributions for the WACS amplitude. 

The perturbative estimate for an IR contribution is given by 

A S (Q) < , (22) 

where summation is over the external lines of the soft subgraph S. The infrared regime corresponds to the Feynman 
mechanism. However, we should keep in mind that the perturbative estimate implies a pointlike coupling of three 
quarks to the proton field while in real life the proton wave function is much softer. In particular, the perturbative 
estimate of the IR regime for the proton form factor gives F(Q 2 ) < Q~ , allowing for 1/Q behavior in principle. 
To get such an asymptotic behavior from our ND models, we should assume that f(x) ~ 1 — x for x close to 1. 
More realistic functions dictate a faster decrease of F(Q 2 ) in the asymptotic Q — > oo limit. No wonder: the IR 
regime is essentially nonperturbative and the Q-dependence of the soft contributions should be better taken from a 
reasonable model rather than from perturbation theory. Again, since accessible Q's are far from being asymptotic, 
the "observed" power behavior of the soft contribution in this region may strongly differ from the asymptotic 
powers given by the Feynman mechanism. 




a) b) c) d) 



FIG. 5. Terms having O(s ) behavior for large s. 

The simplest contributions for the WACS amplitude are given by the s- and u-channel handbag diagrams Fig.|^a,b. 
They correspond to a combined SD-IR regime of Fig.[|a: the dependence on s (or u) is killed by contracting into 
point the quark line connecting the photon vertices while the ^-dependence is killed by removal of a soft subgraph S. 
The SD-regime in this case gives Ajj(s) ~ s° behavior. The nonperturbative part is given by the proton nonforward 
DD's which determine the i-dependence of the total contribution. Another O(s ) configuration is shown in Fig.||c. 
In this case, a hard gluon propagator is convoluted with the proton and photon DD's. Similarly to the usual 
photon structure functions, the photon DD's can be divided into the perturbative and the nonperturbative part. 
The latter corresponds to hadronic component of the real photon while the first one to a direct pointlike quark- 
photon coupling. It can be treated as a part of the one-loop correction to the handbag diagram (see Fig.||d) and is 
accompanied by the a s /ir suppression factor. The hadronic component of the photon DD's has also an extra form 
factor type suppression ~ m 2 /t. 

Just like in the form factor case, the contribution dominating in the formal asymptotic limit s, \t\, \u\ — > oo, 
is given by diagrams corresponding to the pure SD regime, see FigJ^d. The hard subraph H involves two hard 
gluon exchanges which results in suppression by a factor (a s /7r) 2 ~ 1/100 absent in the handbag term. The 
total contribution of all hard configurations was calculated by Farrar and Zhang and then recalculated by 
Kronfeld and Nizic |5(|. Again, a sufficiently large contribution is only obtained if one uses humpy DA's and 1/fc 2 
propagators with no finite-size effects included. Even with such propagators, the WACS amplitude calculated with 
the asymptotic DA is negligibly small pjj] compared to existing data. Our arguments concerning the reliability of 
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CZ enhancements for form factors can be applied to the wide-angle Compton scattering without any changes. For 
these reasons, we ignore hard contributions to the WACS amplitude as negligibly small. 




a) b) c) 



FIG. 6. Configurations involving double and single gluon exchange. 

Another type of configurations containing hard gluon exchange corresponds to the version of the combined SD-IR 
regime shown in Fig.|]c. In particular, they include diagrams like Fig.gb and also diagrams with photons coupled 
to different quarks ("cat's ears", Fig.^Jc). Such contributions have both higher order and higher twist. This brings 
in the a s /ir factor and an extra 1/s suppression. The latter is partially compensated by a slower fall-off of the 
four-quark DD's with t since only one valence quark should change its momentum. 

VI. MODEL FOR WIDE-ANGLE COMPTON SCATTERING AMPLITUDE 



In this paper, we neglect all the suppressed terms and deal only with the handbag contributions Fig||a,b in 
which the highly virtual quark propagator connecting the photon vertices is convoluted with nonforward proton 
DD's parametrizing the overlap of soft wave functions. Since the basic scale 4A 2 characterizing the i-dependence 
of DD's in our model is 2.8 GeV 2 , while existing data are all at momentum transfers t below 5 GeV 2 , we deal with 
the region where the asymptotic estimate (Feynman mechanism) for the overlap contribution is not working yet. 
In the coordinate representation, the sum of two handbag contributions to the Compton amplitude can be written 
as 

M^(p,p';q,q') = f e -*«*V |^ (z/2) 7 ^ c (z) 7 > (-2/2) + ^ a (-z/2) 7 ^ c (-z) 7 ^ a (z/2))|p) d 4 z (23) 

a ^ 

where Q = (q+q')/2 and S c (z) — iz/2ir 2 (z 2 ) 2 is the hard quark propagator (throughout, we use the "hat" notation 
z = z a ^ a ). The summation over the twist-0 longitudinal gluons adds the usual gauge link between the tp,tp fields 
which we do not write down explicitly (gauge link disappears, e.g., in the Fock-Schwinger gauge z a A a (z) = 0). 
Because of the symmetry of the problem, it is convenient to use P = (p + p')/2 (cf. Q) andr=p — p' as the basic 
momenta. Applying the Fiertz transformation and introducing the double distributions by 



{p'\i> a {-z/2)z^ a (z/2)\p) =u{p')zu{p) / dx 



V2 r 



5/2 



-2/2 



< kz *>F a (x, y; t) - e t{kz) F a {x, y; t) 



dy 



+ - u(p')(zf — fz)u(p) I dx 



-2/2 



e -i{k^K a ( x , y; t) - e t{kz) K a {x, y; t)] dy + 0{z 2 ) terms 



(24) 



(we use here the shorthand notation k = xP + yr) and similarly for the parton helicity sensitive operators 



(P' I -0a(-2/2)z 7 5^a(z/2) \ p) = u(p')zj 5 u(p) 



dx 



2/2 



2/2 



e- l{kz) G a {x, y; t) + e t(kz) G a (x, y; t) 



dy 



( rz ) — t i\ f \ 
H u(P nnHP) 



dx 



2/2 



2/2 



e- t(kz) P a (x, y; t) + e l[kz) P a {x, y; t)] dy + 0(z 2 ) terms . 



(25) 



we arrive at a leading-twist QCD parton picture with tilded DD's serving as functions describing long-distance 
dynamics. The new DD's F a (x, y; t), etc., are related to the original DD's F a (x, y; t) discussed in Section II by the 
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shift y — y + x/2. Integrating F(x, y; t) over y one obtains the same nonforward densities T{x\ t). The hard quark 
propagators for the s and u channel handbag diagrams in this picture look like 

xP + yf + Q _ xP + yf + Q xP + yf - Q _ xP + yf - Q 

(xP + yr + Q) 2 xs — (x 2 /4 — y 2 )t + x 2 m 2 (xP + yr — Q) 2 xu — (x 2 /4 — y 2 )t + x 2 m 2 ' 

respectively. We denote s = 2(pq) = s — m 2 and u = —2(pq') — u — m 2 . Since DD's are even functions of y |58| , the 
yf terms in the numerators can be dropped. It is legitimate to keep 0(m 2 ) and 0(t) terms in the denominators: 
the dependence of hard propagators on target parameters m 2 and t can be calculated exactly because of the effect 
analogous to the ^-scaling in DIS || (see also |3(|). Note that the t-correction to hard propagators disappears in 
the large-£ limit dominated by the x ~ 1 integration. The ^-corrections are the largest for y = 0. At this value and 
for x = 1/2 and t = u (cm angle of 90°), the i-term in the denominator of the most important second propagator 
is only 1/8 of the u term. This ratio increases to 1/3 for x = 1/3. However, at nonzero y-values, the i-corrections 
are smaller. Hence, the ^-corrections in the denominators of hard propagators can produce 10% — 20% effects and 
should be included in a complete analysis. In the present paper, we will consider a simplified approximation in 
which these terms are neglected and hard propagators are given by y-independent expressions (xP + Q)/xs and 
(xP + Q)l xu. As a result, the y-integration acts only on the DD's F(x,y;t) and converts them into nonforward 
densities !F(x,t). The latter would appear then through two types of integrals 

/ F a (x,t)dx = F?{t) and f F a (x, t) — = R^t), (27) 
Jo Jo x 

and similarly for JC,Q,V. The functions F*(t) are the flavor components of the usual Fi(t) form factor while 
Ri(t) are the flavor components of a new form factor specific to the wide-angle Compton scattering. In the formal 
asymptotic limit \t\ — » oo, the x integrals for F" (t) and i?"(i) are both dominated in our model by the x ~ 1 
region: the large-t behavior of these functions is governed by the Feynman mechanism and their ratio tends to 
1 as \t\ increases (see Fig.^a). However, due to large value of the effective scale 4A 2 = 2.8 GeV 2 , the accessible 
momentum transfers t < 5 GeV 2 are very far from being asymptotic. 

In Fig.|b we plot T u {x]t) and T u {x\ t)/x at t = -2.5 GeV 2 . It is clear that the relevant integrals are dominated 
by rather small x values x < 0.4 which results in a strong enhancement of R\{t) compared to F±{t) for |i| < 5 GeV 2 . 
Note also that the (p' \ ■ ■ ■ xP . . . \p) matrix elements can produce only t as a large variable while (p 1 \ . . . Q . . . \p) 
gives s. As a result, the enhanced form factors i?"(f) are accompanied by extra s/t enhancement factors compared 
to the Fi(i) terms. In the cross section, these enhancements are squared, i.e., the contributions due to the non- 
enhanced form factors Ff(t) are always accompanied by t 2 /s 2 factors which are smaller than 1/4 for cm angles 
below 90°. Because of double suppression, we neglect F" (t) terms in the present simplified approach. 




FIG. 7. a) Ratio R?(t)/F?(t); b) Functions F*(x;t) (solid line) and J ru (x;t)/x (dashed line) at t = -2.5 GeV 2 . 

The contribution due to the K. functions appears through the flavor components F% (t) of the F-i (t) form factor 
and their enhanced analogues i£§(^)- The major part of contributions due to the /C-type ND's appears in the 
combination R\{t) - {t/ Am 2 )Rl{t). Experimentally, F 2 (t)/F 1 (t) « lGeV 2 /|£|. Since R 2 /F 2 - Ri/F 1 - l/(x), 
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i?2(t) is similarly suppressed compared to Ri{t), and we neglect contributions due to the R%{t) form factors. We 
also neglect here the terms with another spin- flip distribution V related to the pseudoscalar form factor Gp{t) 
which is dominated by the i-channel pion exchange. Our calculations show that the contribution due to the parton 
helicity sensitive densities Q a is suppressed by the factor t 2 /2s 2 compared to that due to the T a densities. This 
factor only reaches 1/8 for the cm angle of 90°, and hence the Q a contributions are not very significant numerically. 
For simplicity, we approximate Q a (x,t) by T a (x,t). After these approximations, the WACS cross section is given 
by the product 



da 



2ira 2 



(pq) (pq') 



(pq') (pq) 



Rl(t), 



(28) 



of the Klcin-Nishina cross section (in which wc dropped 0(m 2 ) and 0(m 4 ) terms) and the square of the R\(t) form 
factor 



Ri(t) 



E- 2 



2 a [RUt) + RUt)] ■ 



(29) 



In our model, R\ (t) is given by 



Ri(t) 



f 

Jo 



el fl a \x) + e 2 fT\x) + 2(e 2 u + e 2 + e 2 s ) (x) 



C xt/Ax\ 2 dx 
X 



(30) 



We included here the sea distributions assuming that they are all equal f sea (x) — f^ s (x) — f a ,d,s( x ) an d using 
a simplified parametrization 



f sea (x) =0.5x-°- 75 (l-x) 7 



(31) 



which accurately reproduces the GRV formula for Q 2 ~ 1 GeV 2 . Due to suppression of the small- a; region by the 
exponential exp[xi/4xA 2 ], the sea quark contribution is rather small (~ 10%) even for — t ~ 1 GeV 2 and is invisible 
for -t > 3 GeV 2 . 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 

-t (GeV 2 ) 



FIG. 8. WACS cross section versus t: comparison of results based on Eq.(28) with experimental data. 
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Comparison with existing data fl59| is shown in Fig.||. Our curves follow the data pattern but are systematically 
lower by a factor of 2, with disagreement becoming more pronounced as the scattering angle increases. Since we 
neglected several terms each capable of producing up to a 20% correction in the amplitude, we consider the agree- 
ment between our curves and the data as encouraging. The most important corrections which should be included 
in a more detailed investigation are the ^-corrections in the denominators of hard propagators and contributions 
due to the "non-leading" K,,Q,V nonforward densities. The latter, as noted above, are usually accompanied by 
t/s and t/u factors, i.e., their contribution becomes more significant at larger angles. The i-correction in the most 
important hard propagator term l/[xu — (x 2 /4 — y 2 )t + x 2 ra^\ also enhances the amplitude at large angles. 

The angular dependence of our results for the combination s 6 (da/dt) is shown on Fig.^. All the curves for initial 
photon ehergies 2,3,4,5 and 6 GeV intersect each other at 8 cm ~ 60°. This is in good agreement with experimental 
data of ref. |59| where the differential cross section at fixed cm angles was fitted by powers of s: da/dt ~ s^ n ^ 
with n cx P(60°) = 5.9 ± 0.3. Our curves correspond to n soft (60°) » 6.1 and n soft (90°) « 6.7 which also agrees with 
the experimental result n cx P(90°) = 7.1 ± 0.4. 




-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 

-Cos(O) 



FIG. 9. Angular dependence of the combination s 6 (da/dt). 

This can be compared with the scaling behavior of the asymptotic hard contribution: modulo logarithms con- 
tained in the a s factors, they have a universal angle-independent power n hard (0) = 6. For 9 cm — 105°, the 
experimental result based on just two data points is n exp (105°) = 6.2 ±1.4, while our model gives n soft (105°) w 7.0. 
Clearly, better data are needed to draw any conclusions here. 
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FIG. 10. s-dependence of the combination s 6 da/dt for = 60° (dotted line), 6 = 90° (dashed line) and 8 = 105° (solid 
line) . 



VII. SUMMARY AND CONCLUSIONS 



In this paper, we introduced nonforward parton densities J-(x; t) which are the simplest hybrids of the usual 
parton densities and hadronic form factors. We proposed a simple model for the quark ND's T a {x]t) which, in 
the t — > limit, reproduces the standard parametrizations for the usual parton densities and gives a reasonable 
description of existing data on the Ff (t) form factor in a wide range 1 GeV 2 < — t < 10 GeV 2 of momentum 
transfer. The crucial observation is that though our model includes only the soft contribution, the form factor 
is dominated at accessible energies by rather small momentum fractions x ~ 0.5 and asymptotic estimates for 
soft contributions (corresponding to Feynman mechanism, i.e., dominance of the x ~ 1 region) are not working 
yet. We gave arguments that the wide-angle Compton scattering amplitude in the same t region is dominated 
by two handbag diagrams. We found also that the largest term contains the same ND's !F(x;t) which determine 
the behavior of Ff(t). However, due to the extra 1/x factor and small value of (x), the WACS amplitude gets a 
strong enhancement bringing our predictions close to existing experimental data. Still, there remains a systematic 
difference by a factor of 2 between our results and existing data. 

On the experimental side, data of higher quality are needed. They are expected from a future experiment at 
Jefferson Lab |30|], in which better statistical accuracy is aimed and several new ideas will be used to control the 
systematic errors. 

On the theoretical side, a more detailed approach is needed which would take into account all nonforward 
densities. A more complete analysis should also include calculable t- and m 2 dependence of the hard quark 
propagators and terms which are not enhanced by the 1 / (x) factors. It should be emphasized that keeping the 
i-terms in the denominators of hard propagators requires a major change in the whole approach: it would be no 
longer possible to get a simplified description in terms of the nonforward densities J-(x;t): one should deal then 
with double distributions F (x, y\ t) in all their complexity and construct a model for their profile in the y direction. 
This observation also demonstrates that the double distributions F{x 1 y; t) are the primary objects for analysing 
nonforward matrix elements of light cone operators. They are more fundamental than their reductions such as 
nonforward, off-forward, etc. distributions which work only when the hard part of the relevant amplitude depends 
on a particular linear combination x + yC, of its two arguments x and y. A more detailed discussion of double 
distributions will be given in a forthcoming publication pl[ . 
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